Long-time behavior of many-particle quantum decay 
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While exponential decay is ubiquitous in Nature, deviations at both short and long times are 
dictated by quantum mechanics. Non-exponential decay is known to arise due to the possibility of 
reconstructing the initial state from the decaying products. We discuss the quantum decay dynamics 
by tunneling of a many-particle system, characterizing the long-time non-exponential behavior of 
the non-escape and survival probabilities. The effects of contact interactions and quantum statistics 
are described. It is found that whereas for non-interacting bosons the long-time decay follows a 
power-law with an exponent linear in the number of particles N, the exponent becomes quadratic 
in N in the fermionic case. The same results apply to strongly interacting many-body systems 
related by the generalized Bose-Fermi duality. The faster fermionic decay can be traced back to 
the effective hard-core interactions between particles, which are as well the decaying products, and 
exhibit spatial anti-bunching which hinders the reconstruction of the initial unstable state. The 
results are illustrated with a paradigmatic model of quantum decay from a trap allowing leaks by 
tunneling, whose dynamics is described exactly by means of an expansion in resonant states. 
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I. INTRODUCTION 

The exponential decay law of unstable systems is found 
across all fields of physics. In the description of nuclear 
stability it was already derived in the early days of quan- 
tum mechanics [1]. However, it is precisely quantum me- 
chanics that imposes deviations of this law during both 
the short and long-time decay [2-.j]. In particular, short- 
time deviations arise from the finite mean energy of the 
initial state and arc associated with the quantum Zeno 
effect, which was experimentally demonstrated in [G] and 
can be exploited for different applications [7]. The long- 
time deviations are exhibited by any physical system, 
described by a hamiltonian h whose spectrum sp(/i) is 
bounded from below. They are generally characterized 
by a power-law decay with a > [■'3, 8]. Elusive for 
about half a century, the observation of this behavior has 
been reported in a recent experiment [9]. Following the 
insight by Ersak, the breakdown of the exponential law 
can be generally attributed to the possible reconstruc- 
tion of the unstable state from the decaying products 
[4]. If this reconstruction is inhibited, the decay dynam- 
ics is governed by an exponential law at all times. This 
argument was recently sharpened in [10] where the inter- 
ference between the reconstructed and non-reconstructed 
state was shown to be responsible for short-time devia- 
tions. Further, it was noticed that the long-time devia- 
tions are indeed due to the state reconstruction in a clas- 
sical, probabilistic sense. In spite of the abundant theo- 
retical work on non-exponential decay, its many-particle 
counterpart is to a good extent unexplored. Indeed, for 
a many-body unstable system the discussion of the de- 
cay dynamics in terms of one-body observables is not 
entirely satisfactory. In the light of the reconstruction ar- 
gument, we might expect correlations between different 
particles arising from quantum statistics and their inter- 
actions to play a crucial role. The short-time deviations 
are under current investigation and will be discussed else- 



where. In this paper, we shall describe the multi-particle 
tunneling decay, paying particular attention to the long- 
time non-exponential behavior and identifying the key 
signatures of quantum statistics and hard-core contact 
interactions. The paper is organized as follows. In sec- 
tion II we introduce the many-particle non-escape and 
survival probabilities, and discuss their computation in 
bosonized and fermionized systems. In section III we 
find their asymptotic behavior in a paradigmatic model 
of tunneling dynamics using the resonant-state expansion 
formalism. The robustness of these results is discussed in 
section IV, and their explanation in terms of state recon- 
struction is provided in section V. We close with a brief 
summary and an outlook. 



II. NON-ESCAPE AND SURVIVAL 
PROBABILITIES OF A MULTI-PARTICLE 
SYSTEM 



Consider a many-body system described by a wave- 
function ^'(xi, . . . , xn; i), symmetric (bosonic) or anti- 
symmetric (fermionic) with respect to permutation of 
particles. For simplicity we consider effectively one- 
dimensional systems relevant to ultracold gases under 
tight transverse confinement where the transverse de- 
grees of freedom are frozen (for this to be the case the 
chemical potential and thermal energy are to be much 
smaller than the transverse vibrational excitation quan- 
tum) [11]. Let us single out a region of interest A, with an 
associate projector (the characteristic function), xa{x). 
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Wc define the N-particle non-escape probabihty as 

N 

-PNit) := (vI/|nXA(x„)|*) 



N 

/ TTda;„|4'(a;i,...,XN;i)r 

JAN 



(1) 



which is the probabihty for the N particles to be found 
simultaneously in the A-region. Indeed, T^n can be ex- 
tracted from the atom number distribution p{n, t) used 
in full counting statistics [12], and defined as the expec- 
tation value 



pin,t) = {6{nA - n)). 



(2) 



where is the density operator in the subspace A and 
n E N. As a matter of fact 



(3) 



where for trapped ultracold gases p{n, t) can be exper- 
imentally measured as in [13]. Hence, Vt^^t) is a truly 
multi-particle observable which should not be confused 
with the non-escape probability defined with respect to 
the density. 



PnW := (^IxaGti)!*) = / dxn{x,t), (4) 

Ja 



where 

n{x,t)= [| d.T„l*(xx2...,a;N;t)l^ (5) 

is the density profile of the cloud, which imposes no con- 
dition on the location of the (N — 1) particles. The be- 
havior of PN(i)i the integrated density profile which is 
a one-body observable, has been the subject of recent 
studies dealing with dynamics of ultracold gases [14, 15]. 

In the description of single-particle quantum decay 
problems of an unstable state ]■(/') (where for N = 1 
the distinction between Vi and Pi becomes superflu- 
ous), the non-escape probability often behaves in 
the same way than the fidelity or survival probability 
Si{t) = |(^(O)|'0(i))|^ The survival probability is re- 
ferred to the overlap between the initial and time-evolved 
state so that its many-particle version simply reads. 



(6) 



and has recently been used to describe Loschmidt echoes 
in one-dimensional interacting Bose gases [Ki]. 



A. Fermionized systems 



ground state is well approximated by a Slater determi- 
nant, 'i>F{xi,...,xt^;t) = ^det^f,^j^(pn{xk;t), where 

(j)n{x,0) is the n— th eigenstate of the initial Hamilto- 
nian, whose time evolution for t > following a quench 
of the trapping potential at t = is denoted by t). 
Using the Leibniz formula it reads ^f{xi, ■ ■ ■ ,xj<s;t) = 

^EpeSN(-l)^nLi'^P(«)(2^",0, where P labels a 
permutation of the symmetric group S'n with N! elements 
and (—1)^ = ±1 is the signature of P. This allows us to 
rewrite the N-particle non escape probability as 

- ^ E(-i)''^'^n(M«)wixAi0Q(„)W) 



P,QGSn 



n=l 
N 



P,R£Sk n=l 



= dctl,^,[{Mt)\xA\Mt))l (7) 

where P,Q,R label different permutations and [17] 

(8) 



{Mt)\XA\Mt)) = / dx0:(.T,t)(/)fc(a;,t). 

J A 

Similarly one can compute the survival probability. 



(9) 



where det^ i^^i[{(l>n{0)\(f>kit))] is the survival amplitude. 
The same expressions in Eq. (7) and (9) describe as 
well hard-core bosons in one dimension, in the so called 
Tonks-Girardeau (TG) regime [is]. Indeed, their wave- 
functions arc related by the Bose-Fermi mapping 



'^tg{xi,- . . ,a;N) = A{xi, xys)-^F{xi, 
through the antisymmetric unit function 

l<j<A:<N 



, a;N), 



(10) 



with e{x) = 1 (-1) if a; > (< 0) and e(0) = . Note 
that this operator is its own inverse and consequently 



and T'^W^^f'^W (11) 



are shared by both dual systems, as well as any other 
multi-particle fermionized system with hard-core inter- 
actions and intermediate statistics, i.e. the whole fam- 
ily of hard-core anyons {"^^hca — A-e'^Fld} for which 
the same mapping holds up to the replacement of A by 



A-6 



<j<k<N ' 



^j) where is the statistical 



parameter [19, 20]. 



Let us consider an atom number state of a spin- 
polarized Fermi gas. Since s-wavc scattering is sup- 
pressed by the Pauli exclusion principle, and p-wave scat- 
tering is generally weak, the Fermi wavcfunction of the 



B. Bosonized systems 

For a non-interacting Bose-gas, the wavefunction 
of the ground state is built as a Hartree product. 
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*B(a;i,...,XN;i) = Yln=i'Pii^ri,t), where (t)i{x,t) de- 
scribes the time-evolution of the single-particle ground 
state of the initial Hamiltonian. It follows that the N- 
particle non-escape probability becomes the N-th power 
of the single-particle non-escape probability Vi (t) , 

N 

^^f^W = {^B\Y[xA\^B)^r,itr, (12) 

where the Vi{t) = J^dx\(f>i{x,t)\'^ is referred to evolu- 
tion of the ground state (f>i{x,t). Similarly, 5^^' = , 
with Si = \{4>i{0)\4>i{t)\'^ . These expressions apply gen- 
erally to bosonized systems, those related by the Bose- 
Fermi duality to non-interacting bosons, like the so- 
called Fermionic Tonks- Girardeau gas and the continu- 
ous family of anyonic extensions [19, 20]. Approaching 
the strictly non-interacting limit in an ultracold atomic 
Bose gas is somehow challenging. The same behavior is 
observed in a weakly interacting gas in the mean-field 
regime, where the Hartree approximation holds, leading 
to 

r^''''^\t)^Vi{tf (13) 

where Vi{t) = J^dx\ip{x,t)\'^ (and analogously for 5n), 
with (p{x^ t) being the solution of the Gross-Pitacvskii 
equation. Nonetheless, an expanding Bosc gas with fi- 
nite interactions eventually acquires a Tonks-Girardeau 
structure [21, 22]. 

III. MODEL 

The study of long-time quantum decay becomes ex- 
tremely challenging or merely intractable with the usual 
numerical propagation methods based of space-discrete 
lattices even with the use of complex absorbing potentials 
in the boundaries of the propagation box [23] . As a result, 
scarce analytic results [3, 24] or specific methods adapted 
to study quantum decay [5] are required. We shall use the 
Resonant-States Expansion (RSE) formalism developed 
by Garci'a-Calderon and coworkers [25, 27], in combina- 
tion with asymptotically exact expansions. We choose 
the one dimensional analogue of the Winter model, one 
of the paradigmatic models of quantum decay by tunnel- 
ing [3] . It consists of an initial box-like potential located 
in the interval [0, a] with infinite walls for t < and whose 
right wall, located at a; = a, is weakened for i > to a 
delta function potential v{x > 0) = i]'5{x~a) of strength 
r] = IrriT]' I'f? > 0, while the left wall v{x < 0) remains 
infinite for all t, restricting the dynamics to x > 0. Its 
piecewise definition allows for a nearly analytical treat- 
ment which facilitates the study of asymptotics both a 
long and short times. Let A := [0,a) be the region of 
interest. We shall study the time evolution of the n-th 
cigenstatc of the initial box-like trap, with the general 

form 4>n{x,t = 0) = J '^sin(knx)x\{x), with fc„ = 



n S N. For t > 0, the time evolution of (/)„ reads 

<i)n{x,t)^ / g{x,x';t)<i)n{x' ,Q)dx' , (14) 

J A 

where g{x,x';t) is the retarded Green's function. The 
RSE [25, 27] will allow us to exploit the analytical prop- 
erties of the corresponding outgoing Green's function 
G~^{x, x'; k), after rewriting g{x, x'; t) as 

oo 

g{x,x';t) = ^u,(x)uj(x')e-*'*'=?*/2m 

+ - / G+(a;,a;';/t)e-'^'='*/2™fcdfc. (15) 
TT Jr 

Here, {uj} arc the resonant states with complex eigcn- 

values Ej ~ = Sj — i"fj/'2, lifetime h/^j and real 

energy ej, which obey outgoing boundary conditions at 
X = a and satisfy the Schrodinger equation. The sum 
runs over the set of proper complex poles kj of the scat- 
tering S-matrix lying on the fourth quadrant of the com- 
plex /c-plane with Re (/cj) > ]Im(A:j)]. The outgoing Green 
function G'^{x,x';k) can generally be written explicitly 
in terms of the regular function and the Jost function of 
the scattering problem [2n]. The integral term involving 
it, is responsible for deviations from the exponential law 
at ultrashort or very long times and we shall focus on its 
contribution. The path F can be chosen as the straight 
line on the complex fc-plane Im(fc) = — Re(fc), passing 
through the origin k = [25]. This expansion holds as 
long as x,x' S A, a limitation that can be overcome 
following [2(>], should that be necessary. Furthermore, in 
the RSE of g(x, x'; t) the sum over resonant states decays 
exponentially with time. Rewriting the countor integral 
as 

- / G+(a:,a:';fc)e-^''''-'*/2™fcdfc 

^ (l + ») (2,s+l) / m X s+k 

^(s-1)! \ht) 

(16) 

where G^'j, = 9^G+(a;, a;' ; fc)|A:=o, we are ready to discuss 
the long-time behavior of both V^(t) and S^(t). The 
long-time asymptotics of the single-particle wavcfunc- 

tion is found to be ^„(z, t) ^ ^^S^ ^^^if^:^ {^f'"' + 

0(t-5/2). leading order, we find the matrix elements 

so that for non-interacting bosons 
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Hence, the exponent of the power-law is hnear in the 
particle number N. The late exponential quantum decay 
preceding the regime in Eq. (19) is governed by the reso- 
nance with longest life time 71 , leading to a multiparticle 
decay rate 



Tn = -d\n'Pn{t)/dt ~ N71. 
Deviations set in during the crossover 
exp(-rNi) « 



(20) 



(21) 



characterized by a non-monotonic behavior, and hence 
the onset of the (5n) asymptotics is the same as in 
the single-particle dynamics, and can be characterized 
by the single resonant criterion, which demands a small 
ratio i? = £1/71 for deviations to occur [2!)]. 

In the case of fermions (N > 1) the first order of the 
expansion vanishes due to symmetry in Eqs. (7) and (9) 
to 0{t~^^)^ and higher order terms are to be taken into 
account. The exact expression of G'^{x,x' \k) in the ID 
analogue of the Winter model can be constructed [30]. 
One can use it to compute the matrix elements enter- 
ing in the definition of 5n and Pn and find the leading 
contribution for different N, which takes the form [31] 



1 



ht 
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{F) 



37r3 (1 + ?7a)4 

3 1 
1757rio [l + 'qaf \ ht 

1024 1 
60153806797r2i (1 + 770)12 ^ ht 

2^ N(2N+1) 



21 



(1 + jya) 



4N 



ht 



(22) 



where the scaling in the last equation has been verified 
for finite values of N by explicit expansion of V^. The 
survival probability exhibits the same long-time asymp- 
totic law, the first few terms being 



Si 



5f) 



1 



ht 



ttS (1 + r;a)4 
729 1 fma^ 

8000000 1 
53144l7r39 + \ ht 

I /^,„2\ N(2N+1) 



(1 + 77a)4N 



ht 



(23) 



The upshot is that for spin-polarized fermions, and 
fermionized systems such as the Tonks- Girardeau gas, 
the exponent of the long-time power-law decay becomes 
quadratic (instead of linear) in the number of particles 
N. 



IV. STABILITY OF THE SCALING 

The leading contribution of every term in the Leibniz 
expansion of the determinant formula of the non-escape 
probability goes as 1/t^, and as l/t^/^ in the case of the 
survival amplitude. The appearance of any power-law 
other than that in Eq. (19) results from the cancela- 
tion of this leading term by symmetry, and the contri- 
bution of higher-order terms. The question arise as to 
whether the symmetry of the system is responsible for 
the multi-particle scaling. The asymptotic behavior is 
robust against a finite shift of the potential barrier to a 
position d 7^ a, d < 00, up to the numerical factors and 
the role played by the initial width of the cloud a. For 
bosons Eq. (19) holds with 



1 



( 



mad 



(l-f77rf)4N \^ ht 
while in the fermionic case 



3N 



N 



1 



(1 + ?7d)4N \ ht J 



(24) 



(25) 



The limit of a vanishing barrier is non-singular and can 
be solved analytically, see Appendix A. The case N = 2 
was recently considered in [33], see also [34]. 



V. DISCUSSION 

One might be tempted to conclude that the faster long- 
time decay of fermionized systems in comparison with 
bosonized gases results from the higher mean energy of 
the former. To appreciate that this is not the case, con- 
sider an excited state of N bosons degenerate with the 
ground state of the non-interacting spin-polarized Fermi 
gas, with one single atom in each single-particle eigen- 



state, ^ft^^; 



(xi,.. 

N 



/N! 



PCi„,fe=i0-i(a;fe;i) 



)p(„-) (a;„, t), EB standing for excited 



^ SpeSN riri^i 'PP{r^ 
Bose gas. The corresponding non-escape and survival 
probabilities can be found to be 



peTl,^,[{Mt)\XA\Mt))l (26) 
\peTl,^,[{MO)\Mm\' ■ (27) 



The density profile of fermions and EB becomes indistin- 
guishable 



(a;,t)=n(^^)(a;) = ^|0fc(x,t)l2, (28) 



k=l 



whence it follows that these systems share the same one- 
body non-escape probability derived from the integrated 
density profile 



(29) 



5 



Nonetheless, 



On 



at variance with 



Vlf^ oc Slf^ oc l/tN(2N+i) 



(30) 



(31) 



Indeed, the same scahng in Eqs. (19) and (30) holds for 
thermal states in the canonical ensemble, showing the 
independence on the mean energy of the initial state. 
Moreover, these long-time power-laws are not specific of 
the ID Winter model, but actually holds for all potentials 
in half-space, where the RSE formalism can be applied 
[25, 27]. 

We further note that from the invariance of T^n and iSn 
under Bose-Fermi mappings, it is actually inappropriate 
to attribute the different scaling to the symmetrization 
imposed by quantum statistics. The different asymptotic 
behavior is indeed a consequence of the hard-core inter- 
actions, which lead to fermionization in the sense of the 
Bose-Fermi duality [19], not to be confused with its dy- 
namical counterpart [20, 35]. For polarized fermions, the 
effective hard-core constraint follows from the Pauli ex- 
clusion principle. Either way, these systems exhibits spa- 
tial antibunching which hinders the reconstruction of the 
initial state in comparison with the case of bosonized sys- 
tems, those dual to non-interacting bosons, and free-of 
hard-core interactions. An exponential decay law holds 
for all times in the absence of state reconstruction [4], 
while a slowly decaying asymptotic power-law is favored 
in those bosonized systems which tend to spatial bunch- 
ing, hence facilitating the reconstruction of the initial 
state localized in the trap. Hard-core interactions lead 
to a faster power-law decay as a result of spatial anti- 
bunching. 

In view of the interaction effects on state reconstruc- 
tion just discussed, and the results by Taniguchi and 
Sawada (restricted to free decay in half space of N = 2 
fermions with both attractive and repulsive Coulomb in- 
teractions) [33] , we notice that the stages of quantum de- 
cay can be further modified whenever the inter-particle 
potential is of finite range. In the attractive case, the 
asymptotic power-law can be slowed down with respect 
to that of bosonized systems, while in the repulsive case 
long-time deviations can be suppressed, extending the 
regime of validity of the exponential decay law. 



VI. CONCLUSIONS AND OUTLOOK 



also shown that this behavior is not related to the higher 
mean energy of the initial state, whose effect is restricted 
to the short-time dynamics. As an outlook, we point out 
that the use of the exact Bose-Fermi mapping for one- 
dimensional Bose gases with finite-interactions [22, 3(i], 
might lead to novel asymptotic regimes of tunneling de- 
cay, not captured by mean-field theories. 

Note- After the submission of this work, Ref. [37] has 
discussed the quantum decay in the same potential model 
in the two-particle case. 

It is a pleasure to acknowledge discussions with J. G. 
Muga, M. Pons, A. Ruschhaupt and D. Sokolovski, and to 
thank the Centre of Excellence for Quantum Engineering 
and Space- Time Research (QUEST). 



Appendix A: Free decay in the iialf-axis 

For a vanishing barrier 77 — > 0. the dynamics is free 
in the positive semiaxis and a fully analytical solution 
is available, writing the propagator by the method of 
images in terms of the free one [27] 



go{x,x';t) 



m 
2iTiht 



exp 



l\2 



im{x — x') 
2th 



as 



g{x, x'\ t) = go{x, x'; t) - go(-x, x'; t). 



(Al) 



(A2) 



The free expansion in the whole space, under go(x,x';t) 
was discussed in [20, 27, 32]. Using it in the superposition 
principle Eq. (14) and Eq. (A2), introducing r = ht/m, 
one finds 

4>n{x,t) = —= a/3[M{f3x - a,akn,T) 

+M{l3x,akn,T)]. (A3) 
We have used the definition of the Moshinsky function 



M{x, k,T) 



e 2t 



w{-z), 



where 



(A4) 



(A5) 



and w is the Faddeyeva function w{z) = e ^ erfc(— zz) 



Summing up, we have analyzed the many-particle 
quantum decay from a leaking trap by tunneling, char- 
acterizing the fundamental deviations from the exponen- 
tial decay law. The long-time behavior, governed by a 
power-law in time, exhibits a signature of quantum statis- 
tics and short-range contact interactions which drasti- 
cally modifies the exponent of the power-law as a func- 
tion of the particle number. Hard-core interactions lead 
to spatial antibunching which hinders the initial state 
reconstruction, and ultimately induces a faster power- 
law decay than in those systems free of them. We have 



Either using the asymptotic form of the Moshinsky 
function [27] or the expansion 

g{x,x';t) = -i^ (™)'/\:,' + 0(i-5/2)^ (^6) 

one can recover for bosons the result in Eq. (19). For 
fermionized systems one needs to keep higher orders in 
the expansion, and take into account the slowest power 
law, which leads to the scaling in Eq. (22). 
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